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Common gastropods such as snails crawl on a solid substrate by propagating muscular waves of 
shear stress on a viscoelastic mucus. Producing the mucus accounts for the largest component in the 
gastropod's energy budget, more than twenty times the amount of mechanical work used in crawling. 
Using a simple mechanical model, we show that the shear-thinning properties of the mucus favor 
a decrease in the amount of mucus necessary for crawling, thereby decreasing the overall energetic 
cost of locomotion. 



I. INTRODUCTION 

Common gastropods crawl on land by propagating waves of shear stress, driven by alternatingregions of muscular 
contraction and expansion, on top of a thin film of viscoelastic mucus released from their foot pIEl Q| The only 
interactions between the organism and the substrate are through this thin mucus layer, hence locomotion is made 
possible solely via the fluid dynamics within this thin film. The nonlinear rheological properties of the mucus are 
responsible for the remarkable ability of these mollusks to walk on solid ground without detaching the foot from the 
substrate Mucus production accounts for about a third of the total energy budget of the animal, and is an order 

of magnitude larger than the mechanical work required for locomotion, making slug crawling the most energetically 
expensive mode of locomotion known among vertebrates and invertebrates |7l l8l l9Ul0|. In this paper, we address the 
relationship between the cost of locomotion and the mechanical properties of the mucus. 

Locomotion strategies employed by gastropods have been of interest to the scientific community for more than 
a hundred years, starting with the works of Simroth 11], Vies an d Parker 01; the reader is referred to Refs. 
|lU2il.ll l Ul llll r J . ll(| for reviews. Gastropods possess a single flexible foot enhanced with arrays of cilia. In many cases, 
the motion of the cilia is responsible for movement of the mucus underneath the gastropod and the resulting motion 
of the animal. However, many gastropods, including all snails that crawl on hard substrates, use a second mechanism 
to crawl, namely the actuation of muscular pedal waves. Using a series of foot muscles, gastropods are able to stress 
the thin layer of mucus (which is typically on the order of tens of microns thick) with alternating regions of muscular 
contraction and expansion. These regions oscillate in time, leading to traveling waves that shear the mucus, resulting 
in translation of the animal with crawling velocities typically between 1 mm/s and 1 cm/s. The muscular waves are 
said to be direct if they propagate in the same direction as the direction of locomotion, and retrograde otherwise. 
Direct waves are waves of contraction, retrograde waves are waves of extension. Similarly, waves are classified as 
monotaxic when a single wave spans the foot of the animal and ditaxic if if two alternating waves span the foot. In 
the three large classes of gastropods, pulmonates are (mostly) terrestrial crawlers and use direct monotaxic waves. 
Prosobranchs represent the largest class of gastropods and are the most varied in structure; most of them are marine 
gastropods using ditaxic retrograde waves. Finally, opisthobranchs do not crawl, but burrow or swim [17| . 

Lissman was the first to study in detail the kinematics and dynamics of gastropod motion [l8L Il9| . He analyzed 
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the locomotion of three pulmonates, Helix aspera (monotaxic direct), Haliotis tuberculata and Pomatias elegans (both 
ditaxic direct). Jones and Trueman |'2(1 \'2l\ were the first to study retrograde wave locomotion, on Patella vulgata 
(ditaxic retrograte). However, it was not until 1980 with the work of Denny and co-authors that a complete mechanical 
picture for gastropod locomotion emerged. The key feature of gastropod locomotion lies in the mechanical properties 
of the mucus layer. This was first discovered for the Pumonate slug Ariolimax columbianus but is generally valid 

for all gastropods ji3)E3| (see also [3,liHEEIi3)- The mucus is composed of more than 95% water, some dissolved 
salts, and the remaining 3-4% is a high molecular weight glycoprotein (mucin) which, in solution, forms a cross-linked 
gel network. This network is responsible for the elastic component of the mucus. At small shear strains, the mucus 
behaves like an elastic solid, but at large shear strains (on the order of 5), the mucus yields and behaves like a viscous 
liquid, with the shear viscosity dropping by about three orders of magnitude to a value of approximately twenty times 
that of water. After yielding, the mucus locally heals back into a gel network on a time scale of less than 0.1 s [3^. 
Underneath a crawling snail, shear strains vary in time between 0, in the so-called "interwave" regions in which the 
foot does not move relative to the ground, and approximately 50 in the regions of largest relative displacement. It is 
this combination of yield-stress properties with judicious periodic shearing of the mucus that allows parts of the foot 
to adhere to the ground while other portions are displaced, leading to an overall translation of the organism. 

The use of a viscoelastic mucus for locomotion has a second interesting consequence. Measurements of the net 
increase in oxygen consumption due to locomotion for Ariolimax columbianus showed that adhesive crawling is the 
most expensive mode of locomotion known Q, twelve times as costly as running for example. Perhaps surprisingly, 
the largest component of this energy budget does not arise from the mechanical work of the muscles, but is due to the 
chemical cost of mucus production (that is, the chemical cost of creating the glycoprotein based on its sequence of 
amino-acids) . For Ariolimax columbianus, mucus production is about twenty times more costly than the mechanical 
work of locomotion jTj. Since Denny's seminal work, the cost of locomotion has been analyzed for other gastropods, 
including Patella vulgata Buccimum undatum 9] and Haliotis kamtschatkana [2sj and similar results have been 
obtained (see also 10] ) . 

In this paper, using a simple model for the nonlinear rhcology of the mucus, we show that its shear-thinning 
properties allow the gastropod to decrease the aforementioned energetic cost of locomotion by minimizing the amount 
of fluid dispensed. The paper is organized as follows. We introduce our model for the nonlinear rheological properties 
of the mucus in flTTl set up the geometry and notation for the particular mathematical problem to be solved in £111 Al 
and present its asymptotic solution in W Bl We apply it to the typical stress distribution underneath a crawling slug 
in 31111 and show that choosing the mucus to be shear-thinning allows to decrease the amount of required fluid. 



II. MODEL 

Adhesive locomotion, the crawling strategy employed by gastropods, has therefore two distinctive features: (1) The 
nonlinear rheology of the mucus, solid-like at low stresses and liquid-like for larger stresses, and (2) the high cost 
of locomotion due to the internal chemical production of the proteins which compose the mucus. In this paper, we 
address the possible link between these two features. Generically, a gastropod should be able to crawl on any kind of 
non-Newtonian mucus, whether it is shear-thinning (as a dilute polymer solution), shear-thickening (as cornstarch in 
water), or any other rheology characterized by a nonlinear relationship between stress and shear rates. 

The mucus of real snails possesses a yield stress, as demonstrated by the experimental work of Denny [5j,|6|, and 
is therefore abruptly shear-thinning with shear viscosities dropping by more than three orders of magnitude around 
the yield point [29j. Here we address the question: Could it be that this particular rheology has been tuned to allow 
locomotion at the lowest energetic cost to the animal? To understand the evolutionary process that could lead to such 
a mucus, consider an organism mechanically similar to a snail that secretes a simple lubricating fluid (e.g. water). 
This mucus becomes more complex as its chemical composition evolves. Hence, as its rheology is perturbed away 
from Newtonian, it is necessary to elucidate which changes to the rheology are beneficial to the organism. This is 
the problem we propose to address, perhaps shedding some light on the importance of material properties in an 
evolutionary context. 
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FIG. 1: Schematic representation of the model used to represent a slightly non-Newtonian mucus in this paper. The straight 
line represent a Newtonian mucus; when the stress-shear rate curve is above the Newtonian one (e > 0) the mucus is shear- 
thickening, and when it is below (e < 0) the mucus is shear-thinning. 

Guided by this evolutionary picture, we make the following modeling assumption: Let us consider the mucus of this 
primitive snail-like organism as slightly non-Newtonian, as illustrated in Fig.^ with a quadratic relationship between 
shear rate, 7 (= du/dy), and shear stress, a, given by 



As a first approach to the problem, we model therefore the mucus as purely viscous (generalized Newtonian) and 
neglect its elastic properties. Our empirical quadratic model (Eq. [Q, a Cross model of order one, is equivalent to 
a Taylor expansion at small stress of the inverse of the shear viscosity, r\ (= cr/7), as a function of the absolute 
value of the shear stress, and is arguably the simplest model that captures the effects of the nonlinear rheology [3t| . 
Here 0% > is a typical shear stress, fi the viscosity in the limit of small stress and e a small parameter, |e| <C 1, 
whose absolute value quantifies the departure from Newtonian behavior and whose sign determines the rheological 
properties (e < corresponds to shear-thinning and e > to shear-thickening) and will be chosen to minimize the 
cost of locomotion. 

As we have seen before, this cost is related to the cost of producing the mucus which differs from the traditional 
measure of energy expenditure associated with mechanical work. In real gastropods, the mucus is dispensed from 
membrane-bound vesicles at various location along foot [3o| | . Let us denote by e the cost of fabricating the mucus per 
unit mass, p the mucus density and Q s the flow rate of the mucus flowing underneath the snail and measured in a 
frame moving with the snail. Then the cost of locomotion, W, over a time interval At, is given by 



and minimizing this cost is therefore equivalent to minimizing the required amount of mucus (per time), Q s . 

In the reminder of this paper, we solve the equation of motion for the crawling snail perturbatively in e and find 
which material properties minimize the flow rate of the mucus, thereby minimizing the overall cost of locomotion. 




(1) 



W = epAtQ 



(2) 



A. Equations of Motion 



Consider a snail that crawls in the x-direction on top of a thin liquid film of mucus of constant thickness h (see 
Fig. |2| . At y = h, the gastropod exerts on the mucus, in a frame moving with the snail, a known traveling wave of 
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FIG. 2: Notation for the model snail. Details are given in Al 



shear stress, cr xy (x, h, t) = a(x — ct), where c is the wave velocity as dictated by the muscular contractions. At y = 0, 
the mucus adheres to the substrate and satisfies a no-slip boundary condition. We denote the unknown instantaneous 
crawling speed of the center of mass of the gastropod by V s e x , where e x is the unit vector in the ir-direction. In the 
reference frame moving with the shear wave, the mechanical problem is time-independent. In this frame, equilibrium 
for the incompressible mucus at low Reynolds numbers can be written as 

Vp = V • er, V • u = (3) 

where p denotes the pressure field, u the velocity field and a the deviatoric stress tensor, subject to the boundary 
conditions 

&xy{x- h) = cr(x) and u(x, 0) = —c — V s , (4) 

where u denotes the velocity field. The equations above are to be solved subject to three constraints: (1) The organism 
is force-free; (2) by periodicity, there is no net build-up of pressure over one wavelength of muscular action within the 
mucus layer; (3) in the frame moving at speed V s e x , the instantaneous velocity of the center of mass of the gastropod 
is zero (by definition of V s ) i.e. 

A _ f X dp f x 

cr(a;) dx = 0, / — [x, y) dx = 0, and / [u(x, h) + c] dx = (5) 



io Jo dx j 

where A is the wavelength of the periodic muscular shearing of the mucus. We further simplify the problem by noting 
that ft < A, and therefore we can apply a lubrication approximation. Consequently, the only relevant shear stress is 
<7 xy , which we denote as a for simplicity, and the only relevant velocity component is u — u • e^,. The equations for 
steady state crawling then become 

dp = da dp = Q 

dx dy ' dy 

Using the stress boundary condition, these equations can be integrated once to obtain 

dp _ 

a{x,y) = -r-(y - h) + <t(x)- (7) 
dx 

The shear stress in the mucus, given by Eq. Q, is therefore a linear function of y. We will denote the coordinate 
where the shear stress changes sign by y*{x), that is 

*(x,y) = ^-(y-y*). (8) 
dx 

Note that above y*, the stress has the same sign as a, and below the opposite sign. As will be seen below, the 
leading-order term for y* is constant and equal to h/3, so that the stress does indeed change sign within the mucus. 
Finally, the first two constraints given in Eq. JSJ lead to 

" x dp 

^fy*dx = 0. (9) 
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d ^= 1 -¥(y~y*) + t sen{y *~ y) m ^(y-y*) 2 - (10) 



B. Asymptotic Solution 

It is now possible to solve for the dynamics of the crawling motion. From Eq. (JJJ and JHJ, the equilibrium equations 
to be solved become 

2 

dy /idx va a ' ' ^ Gif \dx ) v{x) 

Using the boundary condition u(x,0) = — c — V s , and the fact that the velocity profile must be continuous at y* , 
Eq. ijTU|) can be integrated to give 

u(x,y) = |;y(y-V) + e0;[(y-?/*) 3 + 2/* 3 ]-c-^, y < y* (11a) 

u(x,y) = ±y(y-2y*)-e^-[(y-y*) 3 -y* 3 ]-c-V s , y > y* (lib) 

where we have defined ^ = g(x) and = S(x) in order to simplify the notation. The mucus velocity at y = h 

relative to the moving gastropod is therefore given by 

u s (x) ee u(x, h)+c=^(h- 2y*) - ep?-[(h - y*) 3 - y* 3 } - V., (12) 
and the flow rate in the frame moving with the gastropod is given by 

9 h2 ( h e 9 2 S r * 4 , , u „,*n4 ,.*3i 



[u{x,y) + c] dy = »— - y*\ - V s h - ^^[z/* 4 + {h- y*T ~ V d *]- (13) 



Finally the average rate of energy dissipated through viscosity (which is approximately equal to the rate of mechanical 
work done by the gastropod to crawl) is given by 

/> A r-h 

£ = / / ajdydx. (14) 
Jo Jo 

The solution can now be expanded in powers of e and solved order- by-order. For details see Appendix [S] Recall that 
our goal is to find the influence of the sign of e on (1) the chemical work required to produce the mucus, which is 
proportional to Q s as defined in Eq. (|13fl . and (2) the mechanical rate work of crawling, as defined in Eq. I|14|) . 
At lowest order, the flow is Newtonian. The constraints are written at this order as 

A /• A t> A 

go dx — and / g^yo da; = 0, / uo(x, h) dx = 0. (15) 
o Jo Jo 



Averaging the zeroth order velocity field (Eq. IA2a|) over one wavelength and applying the constraints leads to a lowest 
order crawling velocity Vq = 0. Thus, as expected, the organism cannot propel itself on a Newtonian fluid. Averaging 
the volume flux fEa. IA2d|) over one wave length, gives Qo = 0, and therefore yo = h/3. Using this result in Eq. © 
we find the pressure gradient and the energy associated with mechanical work are given by 

3ct(x) hX 2 

90 ^^2h = T^' ' ^ ^ 

respectively, where we have denoted (/) = i L f(x)dx = f(Xu)du. 
At order e we have to enforce the constraints 

/•A /*A /*A 

/ 5idx = 0, / [g yx + giy ]dx = 0, / m(x, h)dx = 0. (17) 
Jo Jo Jo 



Averaging the velocity at this order (Eq. IA2b|) leads to a non-zero first order velocity 



Vi = ~-^{a\a\). (18) 
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So in general, the organism can crawl at order e. The motion is retrograde if e{cr|cr|) > 0, and direct otherwise. 
Averaging Q at order e (Eq. IA2cl) then gives the first order mucus flow rate 



79 h 2 
432 jjicr* 



Note that V\Qi < as is expected. Substituting this result back into the first order flux equation (Eq. IA2e|) 
(recall Q s is a constant independent of x) and using both the lowest order solution for yo and Eq l|18|) leads to 
9oVi = — fj|(r|)/(216(T*). Expanding Eq. JSJ at order e gives g\ = 2>goy\/2h, leading to the first order term in 

the pressure gradient 

^ = ~<m)-m' (20) 

Finally, we find the first order correction to the mechanical energy 

9o ^icr* 

Consequently, for a known traveling wave of shear stress applied by the gastropod on the mucus, the least amount of 
mechanical work is done when e > 0, that is when the mucus is shear-thickening. 

Since we find that Q — 0, we have to calculate the second order correction to the flux in order to quantify the 
influence of the sign of e on the chemical cost of locomotion. At this order we have to enforce the constraints 

A p A r A 

g 2 dx = 0, / [g y 2 + giV® + giyi]dx = 0, / u 2 (x, h)dx = 0. (22) 



Again, averaging the velocity at order e 2 (Eq. IA2cl) we find the second order correction to the crawling velocity 

5 h 



V 2 = J^^(\°\)(°\°\}-(° 3 )1 (23) 
486 /laf 

and averaging Q at order e 2 (Eq. IA2f|) gives the second order correction to the flux 



III. RELEVANCE TO GASTROPOD CRAWLING: IMPACT OF RHEOLOGY ON LOCOMOTION COST 

We are now equipped to address the main question raised in this paper, as presented in iflTl namelv. have the material 
properties of pedal mucus been tuned to optimized crawling efficiencies. In the frame moving with the gastropod, the 
flow rate is given by Q s — eQi + e 2 Q 2 , with errors of order e 3 . This flow rate is, as discussed above, directly related 
to the cost of locomotion associated with the chemical production of the mucus. A change of sign for e modifies the 
sign of the first order term, but does not change the second order term, and therefore the absolute value of the flow 
rate is minimized for a particular choice of the sign of e. Note that we consider the absolute value of the flow rate to 
allow for both direct and retrograde waves; i.e. the direction of motion relative to the direction of the traveling wave 
is irrelevant. The absolute value of the flow rate normalized by the first order term is given by 

l0 -'=l + ^, (25) 



with 



which we estimate below. 
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Let us consider the distribution of shear stress underneath a crawling gastropod, cf(x), and let us zoom in on the 
details of one wavelength A of the foot, as illustrated on the left in Fig. [21 We know that the stress, which averages 
to zero, will be negative in some region of the shear wave and positive in the rest. We can safely assume that it is a 
smooth function of space, and we write for convenience 

a(x) = -a-g(x), for xe[0,aX], (27a) 
a(x) = a + f(x), for x£[aX,X], (27b) 

where / and g are positive functions with maximum value I, a £ [0, 1] represents the fraction of the wavelength which 
has negative stress, and <7_ and <r+ are the maximum absolute values of the regions of negative and positive stresses 
respectively. As indicated in Fig. [21 we have defined the wavelength such that cr(0) = <r(A) = and, given Eq. H27JI . 
the stress is also equal to zero at x — aX. Physically, the stress curve in Fig. [21 can be thought of as the profile of the 
muscular strains in the foot, and therefore the area where the stress a goes from negative to positive (around x = a A) 
is a region of muscular extension and that where a goes from positive to negative (around x — and A) is an area of 
muscular contraction. 

The profile of shear stress exerted by the gastropod on the mucus is a result of the action of its foot muscles. 
The most detailed studies of foot muscles in real snails are due to Jones and coauthors, for Patella vulgata ^2j^ and 
Agriolimax reticulatus 31] . They showed that the distribution of muscle fibers in gastropod varies widely between 
species, as confirmed by subsequent studies of Neritina reclivata and Thais rustica [32|, Helix aspersa [33j, and 
Melampus bidentatus [34| . However, they also showed that (1) the distribution of muscles responsible for the wave 
motion are front-back symmetric in each of the the organisms considered (dorso-ventral muscles for Patella vulgata 
leading to retrograde waves, and longitudinal, transverse and oblique muscles for Agriolimax reticulatus and direct 
waves), and (2) that the same sets of muscles were responsible for both extension and contraction of the foot. Based 
on these observations, we will assume in the remainder of this paper that the action of the muscles result in wave 
profiles in which the regions of muscular extension and those of muscular contraction have similar profiles. In other 
words, up to scaling factors in both the x and y directions, regions of positive and negative stress as sketched in Fig. [21 
will be assumed to have similar shapes. We expect this assumption - the simplest assumption that can be made about 
the stress profile in Fig. [21- although it is restricted to organisms with similar characteristics to Agriolimax reticulatus 
and Patella vulgata, to remain true is other types of gastropods. 

The mathematical formulation of this assumption is given by a relationship between / and g as 

9{x)=f(—x + x\- (28) 



For simplicity we introduce a new positive function, w(u), defined for u € [0, 1] as w(u) — g[aX(l — it)] = f[X(a+ (1 — 
a)u)]. Since the organism as a whole is force free, 

a(x) dx = Xa ~ a - \ h = 0, (29) 

where we denote 

I n = [ w n (u)du > 0. (30) 
Jo 

<r_ = — -a+- (31) 
a 

We can now proceed with the calculation of the various averaged quantities along one wavelength of the foot. To 
calculate the effect of the second order correction to the flux from Eq. I)26|l . we require expressions for {a\a\) : (a 3 ) 
and Rewriting these quantities in terms of cr+, a and /„, we obtain 

{<r\<t\) = -vl p X g 2 (u)du + * 2 + t f 2 (u)du = (2a ~ 1)(1 ~ a) 4/ 2 , (32) 



Since /„ is strictly positive, we require 



nA 
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<a 3 > = 



tX pX 

g 3 (u)du + al / f{u)du 

JaX 



(2a-l)(l-g) 3 

~2 a + l3 ' 



(33) 



(\a\)=a. 



aX pX 

g(u)du + a + / f(u) du = 2(1 — a)cr+Ii. 



lo JaX 

Thus the excess flow rate, as calculated form Eq. Ij26(l . is given by 

~Q~ ~ £ 4266 I cr* J I a 1 '' ' 



(34) 



(35) 



where we have defined f3 = Recall that we are interested in the sign of this term to determine which sign 

of e is beneficial to the cost of locomotion. Hence we need to determine the sign of a 2 — a + (3, as all of the other 
terms on the right hand side of Eq. Il35[l are strictly positive. It is straightforward to show that this term is positive 
if 4/3 > 1 or equivalently, if 2/3 > This inequality can be demonstrated as follows: 



2U 



hh = 



w 2 (x) 



2w(x) — / w(u)du 


1 



da; 



uu (x) 



2w{x) 



w(u)du 



dx 



w 2 (x) 



2w(x) 



w(u)du 



dx 



> - 



w(u)du 

















1/ 


2w{x) — 1 w(u)du 


dx + 


2w(x) — / w(u)du 


dxj 
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w(u)du I > 0. 
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We have therefore shown that a 2 — a + (3 > 0, which implies that the sign of AQ/Q is the same as the sign of e, for 
all possible functions w. Given Eq. I|25|l . this means that the absolute value of the flow rate will be minimized if we 
choose e < 0, i.e. a shear-thinning fluid. As chemical production represents the largest energy cost associated with 
locomotion, this criteria trumps the mechanical work result given above, and a shear-thinning fluid will decrease the 
overall cost of locomotion for the gastropod. 



IV. DISCUSSION 

In this study, we have investigated the possible link between rheological properties of gastropod mucus and the large 
cost of locomotion experienced by these animals. In theory, any type of nonlinear relationship between shear stress 
and shear rate should allow the animal to move, hence there are a number of open questions regarding the possible 
mechanisms and criteria used in selecting particular types of mucus. The approach we take is to use a generalized 
Newtonian model (Fig. ^ and Eq. ^) to quantify how to best perturb the mucus away from a Newtonian fluid. Using 
a series of assumptions which are justified from an experimental standpoint, we have been able to show that a mucus 
which is shear-thinning allows the gastropod to crawl while using the least amount of fluid. As the energetic cost 
of locomotion for these animals is largely dominated by that of mucus production, a shear-thinning material allows 
locomotion at the least expense. This is consistent with experimental observations as all studies which have quantified 
the rheology of real gastropod mucus have found the fluid to be strongly shear-thinning. 

Real mucus is abruptly shear-thinning, so is tempting to suggest that our simple model can provide some rationale 
for the biological series of events that transformed a simple lubricating fluid into an extremely nonlinear material 
through the amplification of its beneficial nonlinear characteristics. Certainly the lowered cost of locomotion is an 
added benefit for the animal regardless of whether or not the energetics played a significant role in the evolutionary 
process. However, our calculations were made under several simplifying assumptions, the most severe of which is the 
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focus on viscous features while ignoring elastic stresses. Consequently, at best, we can conclude that our calculations 
suggest an interesting relation between locomotion cost and the mechanical properties of the fluid. 

The second point of interest is that if, instead of using what we know to be the real measure of locomotion cost for 
snails - the chemical cost associated with mucus production - we use a more traditional measure of locomotion cost - 
the mechanical work of the crawling organism - we find the opposite result. That is, wc find that it is beneficial to use 
a shear-thickening fluid to decrease the amount of work necessary to move. This erroneous conclusion emphasizes the 
importance of biologically relevant measures when addressing optimization and tuning of biological systems, and it 
suggests that the results of our calculations have some biological relevance. Note that it also implies that for synthetic 
crawlers, for which mucus production is not an issue, shear-thickening fluids are energetically advantageous j35| |. 

Our third conclusion concerns the speed of the muscular waves. We have found that crawling speeds and flow rates 
are independent of the speed c of the muscular waves, which might at first appear counterintuitive as it allows for 
locomotion when c = 0. However, we have assumed implicitly in our calculations that the size of the snail foot (and 
therefore its wavelength) remains constant. Thus we have assumed that the muscular extensional strains in the foot 
remain small, as is observed experimentally Q Q. In the frame moving with the snail, the typical stress a is 
sustained by the same sets of muscles during a time t ~ A/c, and the typical local velocity of the foot relative to the 
ground is such that <r ~ fxv/h so that v ~ ha/fi. The maximum longitudinal displacement of elements on the foot is 
therefore I ~ vt ~ ha\//ic, and therefore the typical muscular strain is of order £/X ~ ha/fic. Requiring an upper 
bound on muscular strains offers therefore a lower bound for the wave speed, removing the specious c = paradox. 

Finally, we discuss the difference between direct and retrograde wave crawling. As we have seen, the crawling speed 
is given at leading order, by V s = — 7eh(a\<7\} /36[ia* and, in theory, both direct (V s > 0) and retrograde (V s < 0) 
motion is possible. On the other hand, if the foot of the snail lifts from the substrate as the wave passes - even slightly, 
as in most snails - the contribution of this lifting motion to locomotion always results in a retrograde force due to 
peristaltic effects pifij . The consequence of this result is that, if the snail uses direct crawling yet lifts its foot, the 
contribution due to the non-Newtonian rheology of the mucus needs to be large enough to overcome the retrograde 
contribution due to foot lifting. More precisely, if the foot shape is written as h(x) — ho[l + 6v(x)], where <5 <§; 1 is its 
amplitude and v(x) its dimensionless profile, using the results in Ref. [3^|, we find that the crawling velocity is given, 
at leading order in both e and <5, by 



V. = GcS 2 



7 eh 

— (a\a), 

6b /Lt<7* 



(36) 



and direct wave motion (V s > 0) is possible only if the mucus is sufficiently non-Newtonian, namely 

216 c<5 2 ^ct* 



e(a\a\)<0, \e\ > 



7 h\(a\a\)\ 



(37) 



In the case where direct motion is preferable from a muscular standpoint, the result of Eq. (|37|l suggest an additional 
incentive for the mucus to depart from Newtonian. 

APPENDIX A: ASYMPTOTIC RESULTS TO ORDER e 2 

The solution to Eq. (|10|) is expanded in powers of e as follows: 



' u s (x) 




' u n (x) " 






9n(x) 


< y*{x) 




Vn(x) > 


v s 


n 


v n 






Qn 



(Al) 
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Grouping like orders of e, we find 



u (x) 
u\(x) 
u 2 {x) 



-fjo 



yo - Vo, 



h 


'* 


fh 






fi- 




V2 




-302/1 


ll 




fh 






/' 


92 


V2 


- yo^j 


-3i2/i 



ft 2 /ft 

2// 



Qo = —.9o ( g - 2/o 



ft 2 




fh 






~ 2^t 


3i 


V3 


- yo^j 


-302/1 


/i 2 






- 2/o 1 




~ 2/x 


32 


V3 




-312/1 



+ |^[2/o 3 + (2/o-^) 3 ]-V 1 . 



V 2 

+ (yo - M 3 ] + 3yoyi [y 2 + (yo - hf 
hV , 

9lS 



hVi 



12^(7* 

hV 2 , 



+ (h- y ) 4 - 4fcj/g 



[31 [yo + (ft - yo) 4 - 4fry 3 ] + 2 5o yi [»g + (yo - ft) 3 - 3/i^]] , 



fib 
£1 



1 



^ y 3o[yo + (^-yo) 3 ]da 
1 f A 

3a* Vo 



1 



4/^(7 



[3yi3o[yo - (ft - yo) 2 ] + 23031 [(ft - y ) 3 + yo]] dz 

A 

-3 or ,4 



* JO 



g*S{y 4 Q + (h- yo) 4 ] dx, 



(A2a) 
(A2b) 
(A2c) 

(A2d) 
(A2e) 
(A2f) 

(A2g) 
(A2h) 



which can be solved order- by-order. As we have shown in main body of the paper, it it not necessary to calculate £ 2 
to determine the effect of the sign of e on the mechanical rate of work. 



APPENDIX B: ANALYTIC EXPANSION 

The model presented in the paper is quadratic in the stress-strain rate relationship. We present here an alternate 
model, where the relationship is cubic, the lowest-order model where this rheological relationship is analytic, and is 
written as 

7 = -(l-4V (Bl) 

We use the same notation as in the main part of the paper. The constraints on the motion, as given by Eq. |J5}, 
remain valid here. The leading order solution is the same as in the main part of the paper, and we find 

Vo = 0, (B2) 
Qo = 0, (B3) 

£0 = £V>, (B4) 
ft - |. (B5) 



At next order, we obtain crawling with 
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« - -h^- (B6) 

* - < B7 > 
£ ' = -To^- < B8 > 



The second order solution is then given by 

9 ft 



6400 ixai 
and therefore we have 

AQ = 21e (<x 5 )-(<x 3 )(<x 2 ) 
920<r 2 (a 3 ) 

Now, using the same arguments and notation as in AIIII we find 



^ 2 = 7^H^4^ 3 >^ 2 >-^ 5 >]> ( B1 °) 
1280 /^crj 

21 ^l^ 5 )-^ 3 )^ 2 )], (BH) 



(B12) 



(a 2 ) = a 2 _ g 2 {u) du + a 2 , C f(u) du = (1 "K lh, (B13) 
Jo Ja\ a 

t-B\ 5 r A 5mj 5 / , * rfn . (2a 2 -2a + l)(2a-l)(l-a) 5r 

and therefore 

AQ_ 21 f<7 2 \J 5 ^ 2a 2 -2 Q + l + (a 2 -q)C ^ 

"Q" " £ 920 73 I ? J ' (B15) 

where £ = ^Iz/h- It is straightforward to show that the quadratic numerator in Eq. (|B15|) is positive if and only 
Q < 2, that is, I2I3 < 2/5. This is always true and can be shown in the same fashion as the demonstration that 
hh < 2/3 in STU 

In conclusion, we show that with this model as well, the sign of AQ/Q is the same as the sign of e, and therefore 
the mucus should be chosen to be shear-thinning (e < 0) in order to minimize the amount of mucus production, and 
therefore the overall cost of crawling. Note again that since E\ < 0, this conclusion is opposite to the one we would 
have obtained if mechanical work was our measure of locomotion cost. 
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